Homodyne estimation of Gaussian quantum discord 
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We address the experimental estimation of Gaussian quantum discord for two-mode squeezed thermal state, 
and demonstrate a measurement scheme based on a pair of homodyne detectors assisted by Bayesian analysis 
which provides nearly optimal estimation for small value of discord. Besides, though homodyne detection is 
not optimal for Gaussian discord, the noise ratio to the ultimate quantum limit, as dictacted by the quantum 
Cramer- Rao bound, is limited to about 10 dB. 



Quantum correlations are central resources for quantum 
technology. These tight connections empower the advantages 
shown by the exploitation of quantum coding in applications 
to cryptography, computation and sensing. While at first en- 
tanglement was recognized to be the most peculiar form of 
quantum correlations, novel concepts have been introduced to 
capture either more specific aspects, such as quantum steer- 
ing 1 1 , 2 1, or, to the other end of the spectrum, more general 
occurrences. Quantum discord represents the most success- 
ful attempt to observe quantum features within the current 
picture|[3] HJ: it is related to the fact that quantum informa- 
tion in a bipartite system can not be accessed locally without 
causing an inherent disturbance at a difference with classical 
probability distributions. 

Quantum discord has recently attracted considerable atten- 
tion, due to its possible, yet controversial, usefulness as a re- 
source in mixed-state quantum computing. There exist in fact 
architectures for which an exponential improvement over clas- 
sical resources is obtained |l5]|6l, albeit the entanglement be- 
comes exponentially small QIHl- Discord has then been sug- 
gested as the empowering resource, while following investiga- 
tions contested this interpretation |9|. This debate has stimu- 
lated an intense effort into looking at protocols where discord 
acts a resource: it has been demonstrated that discord does 
play a role in the activation of multipartite entanglement lflOl . 
entanglement generation by measurement lITTI . state merging 
lfT2l . and for complete positivity of evolutions lfT3l[T4l . 

In the experimental test of such proposed connections, the 
comparison of discord with relevant figures of merit is clearly 
connected to the ability of estimating with the best precision 
allowed by a given amount to resources. A key problem is 
then to find optimal strategies, and to understand their fun- 
damental limit introducing proper Cramer-Rao bounds (CRB) 
Il2n 423l. In fact, experimental observation of quantu mdiscord 
has been undertaken either by direct inspection of the density 
matrix ifTsHTOl , or by using a witness ||20| , however, with no 
concern about the optimality of the scheme. 

Optimal estimation of quantum coiTelations has been in- 
vestigated for entanglement ll24l and optimal estimators have 
been experimentally proved to attain the quantum limit for 



different families of qubit states f25l. For the perspective of 
quantum metrology, this is highly nontrivial, since there exists 
no observable directly related to quantum discord. A proper 
estimator is then needed, which might depend on several char- 
acteristic parameters of the quantum state. In such a mul- 
tiparameter problem, finding an optimised detection scheme 
might be hard, and could demand complex experimental ap- 
parata or heavy post-processing of the data. 

In this Letter we demonstrate homodyne estimation of 
Gaussian quantum discord in continuous variable systems 
1261127 1. and compare the achieved level of precision with the 
classical CRB for homodyne detection, and with the quantum 
CRB, which sets the ultimate precision allowed by quantum 
mechanics. We found that although homodyne detection is 
not optimal for Gaussian discord, the noise ratio to the ulti- 
mate quantum limit is limited to about 10 dB. Our findings 
also show how a suitable Bayesian data processing may be 
employed to improve precision, especially in the estimation 
of small values of discord. 

Quantum discord is defined as the difference between two 
quantum analogues of classically equivalent expressions of 
the mutual information in bipartite systems. Its evaluation 
demands an optimization procedure over the set of all mea- 
surements on a given subsystem. For continuous variable sys- 
tem, such minimisation reveals as an extremely complex task; 
however, in the case of Gaussian states, we can conveniently 
restrict the search to Gaussian measurement only ||26l , obtain- 
ing an expression for the Gaussian quantum discord Il26ll27l . 
This sets a lower limit to the discord of the state, and also rep- 
resents an operative figure of merit in those context where, for 
experimental convenience, only Gaussian measurements are 
employed. 

Our investigation is concerned with an important class of 
Gaussian states, i.e. the two-mode squeezed thermal states 
(STS) naturally produced by a non-collinear optical paramet- 
ric amplifier (OPA). If we introduce the two-mode squeez- 
ing operator 6*2(5) — exp ( s{ala\ — aoOi) ) , and the thermal 
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giN„ Nt) - S2{s)iy{Nt) ® i^{Nt)S2{sy, (1) 

and thus can be fully described by the two parameters 
A^s=sinh s and Nt, representing, respectively, the effective 
amount of squeezing photons and thermal photons. In fact, 
spurious effects such as unwanted amplification, result in a 
loss of purity of the squeezed state by thermaUsation, but do 
not affect the Gaussian character of the emission, so the form 
of the density matrix ([TJ provides a fully general description 
of the output of a realistic OPA ||28l . 

For the class of states in Eq. (fill the Gaussian quantum 
discord is given by 

DiN,,Nt) = /i(ki) - 2hiK2) + h{K3) , 

where /i(x) = (a;+l/2) log(a;+l/2)-(x-l/2) log(a;-l/2) 
is the binary entropy and ki == (1 + 2Ns){l + 2Nt), K2 ~ 
(Nt + 1/2), K3 = (1 + N, + Nt)iNt + l/2)/(l + N,+ 
Nt + 2NsNt). We can estimate the discord from Ng and Nt 
as we varied the pump power of our OPA 1291 . For each power 
setting, these two parameters are extracted by the outcome of 
two homodyne detectors, one on each mode, which measure 
pairs of quadratures {Xo,Xi} and {Po,Pi} (Fig.l). From 
these, we can evaluate the four linear combinations 
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where Q^^' and Q'^-* are squeezed quadratures, while Q'-^' 
and Q^^^ are anti-squeezed; in particular Mg measurement 
outcomes are recorded for each one of the four quadratures. 
The corresponding variances, cr^{Qsq) and (J^{Q usq), that can 
be obtained from the experimental data, can be rewritten as 
function of Ns and Nt as follows 

(T'(Qsq/asq) = (1 + 2iV. T 2v/A^,(l + iV,))(l + 2Nt) (3) 




FIG. 1: Conceptual layout of the discord estimation. Our non- 
collinear OPA is based on a nonlinear KNbOa cristal, pumped by 
a frequency-doubled Ti:Sapph laser (repetition rate 800kHz, wave- 
length Ap=425nm, pulse duration 120fs). This produces a STS with 
discord D depending on the pump power, i.e. on their average pho- 
ton number, the two mode are measured by two homodyne detectors 
DHo and DHi. The relative phase between the local oscillators is 
locked by mapping them to different polarisations on the same spa- 
tial mode. In this way, we can record blocks of 20000 values of pairs 
of quadratures for {Xq , Xi } and the same for {Po,Pi}. 



where the probability of obtaining the outcome q, by mea- 
suring the quadrature Q^'^^ is a Gaussian distribution 
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For squeezed quadratures (fc — {1,4}) we substitute c^. — 
cr'^(Qsq), while for anti-squeezed quadrature {k — {2,3}) 
a\ — (J^iQasq)- By means of the Bayes theorem, we obtain 
the a-posteriori probability 

piNs,Nt\X) = ^p{X\Ns,Nt)p„{N,)p„iNt), (6) 
AA= f dN,dNtpiX\N„Nt)po{N,)poiNt). (7) 



The expressions obtained can be then inverted to obtain the 
experimental estimate N"'^^ and N"'^^, along with the relative 
uncertainties a'^iN™'') and cr^(iV|"^). These values can be 
used in the expression for discord f29l to calculate its value 
Z?'"^, and the uncertainty a^{D"'^^). The uncertainties on 
these quantities are then obtained by a Monte Carlo procedure 
|J29l. One can use the same data and refine the estimation 
by using a Bayesian analysis. As described above, each data 
sample corresponds to Mq = 2 ■ 10^ measurement of each 
of the four quadratures. The total sample, thus coiTespond to 
Mt — 4:Mq homodyne outcomes 

;r-lo(l) rM^ o(2) (2) (3) (3) (4) (4). 

^ - i?! ,--.9a/,'9i ,--,1m^i1i '■•'9m,''?! ■'■■^ImJ- 
The overall sample probability can be evaluated as 
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where the Po{Ns) and Po{Nt) are the so-called a-priori prob- 
ability distributions for the two parameters. In our proce- 
dure, we use the results of the inversion estimation to con- 
struct these a-priori distributions. That is, we consider po(-^s) 
and po [Xt) as Gaussian functions with respectively, mean val- 
ues equal to Nf" and N™'^ , and variances equal to a'^{Nf") 
and a'^{Nl^^). Then, we can use the a-posteriori probability 
distribution evaluated as in Eq. (|6]l to obtain an estimate of 
the two parameters and of their variances. In formula, (for 

i = s, t) 
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a\N°^^) = j dN, dNt [N, - A^f ^)^ p{N,,Nt\X) . (9) 

By using the formula of the discord for two-mode STS in 
||29ll and by propagating the errors, we then obtain an estimate 
£)bay ^Qj. j.jjg discord, alongs with its variance <t'^{D^'^^). 



The value of discord depends on both the squeezing and 
thermal photons. Consequently, its estimation is inherently a 
multi-parameter problem, and we have to identify the relevant 
physical parameters to evaluate the correct CRB. In the mul- 
tiparametre scenario, the quantum Fisher information (QFI) 
associated to a vector of parameters A={Aj}o<i<ri is in the 
form of a matrix H. This sets a lower bound on the covariance 
afj^{XiXj) — {Xi) (Xj) after M repetitions on the experiment: 
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In the specific case of our experiment, we can bound the 
uncertainty on the discord D of the states we prepare as: 
CT^(-D) > jj{H^^)dd. While our measurement strategy has 
the advantage of being simple, it is not expected to be opti- 
mal, i.e. to saturate the quantum CRB. In order to assess the 
estimator, i.e. the data processing, we also need to compare it 
to the classical CRB associated to our specific measurement, 
which is analogously described by a classical Fisher informa- 
tion (FI) matrix F. 

In the evaluation of the correct bound, we need a suitable 
parametrisation of the state, so that in the expression ( fTO] ) one 
parameter only actually varies, while the others are kept fixed: 
this can not be the case for the number of thermal and squee- 
ing photons, as both of them change with the pump power. 
Therefore, we need to reshape the QFI matrices for different 
couples of parameters, so to consider those which are more 
directly connected to the experimental conditions. We start 
by considering the first couple Ai = {Ns,Nt}; by using the 
formulas described in ||29l , we obtain 



H'^^^ = diag 



(l + 2A^t)2 



Nsil + iV,)(l + 2Nt + 2Nf) Nt{l + Nt 



(11) 



As explained above, thermal photons appear because of im- 
perfections in the operation of the OR4 and because of loss. 
When the squeezing is not too low, we can reparametrise 
our state by taking in consideration the effective squeez- 
ing strength r, and a parasite amplification with strength jr 
ll30l[3Ti . The overall homodyne detection can be separately 
calibrated, obtaining ?/ = 0.62. Thus we can rewrite the ma- 
trix ( [TT| in terms of the two unknown physical parameters 
•^2 — {^-t} via the expression i?"' ~ Bi2H^^' BJ2, where 
B12 is the transfer matrix for this change of variables I.29J . 
Next, since the physical parameter that changes during our ex- 
periment, resulting in the variation of the amount of discord, 
is the squeezing parameter r (while 7 and 77 can be consid- 
ered to remain constant), we perform the last change of vari- 
able, by considering A3 — {D, 7}. Again the QFI matrix can 
be obtained as H^'^' = i?23-ff S23' ^i^^ the bound on the 
variance for the quantum discord can be easily evaluated as 
described in Eq. ( [T0| . 

We also want to derive the classical CRB for quantum dis- 
cord, that we obtain if we consider as measurement homodyne 
detection of squeezed and anti-squeezed quadratures of a two- 
mode squeezed thermal state. Let us start by considering the 



Fisher information matrix we obtain if we want to estimate the 
two parameters Ai ~ {Ns, Nt} by means of homodyne detec- 
tion on a certain quadrature Q^. Since the state is a Gaussian 
state, the conditional probability distribution of measuring a 
value X, is a Gaussian function, with zero mean, and variance 
(t'^{Q^). By using the formulas in ||29l and evaluaiting some 
Gaussian integrals, one easily obtains the following formula 
for the Fisher matrix elements 



1 
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where A^ — {Ns,Nt}. If one considers to measure the 
squeezed or the anti-squeezed quadratures one obtains the fol- 
lowing FI matrices: 
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If we perform a fixed number of measurements, where half 
of them are done on the squeezed quadratures, and the re- 
maining ones on the anti-squeezed quadratures, the overall 
FI matrix which will give the CRB for the two parameters 
Ai = {Ns, Nt} is obtained as 



f(i) = l(F>'q + F'^'^q) =diag( 



2Ns+2Nf {l + 2Nt)^' 
(13) 



To obtain the CRB for homodyne detection of Gaussian dis- 
cord, we can proceed as we showed for the quantum CRB, 
simply replacing the QFI matrices, with the FI ones. The val- 
ues of the discord obtained using our Bayesian estimation are 
shown in Fig. l2] the points indicate the experimental data, 
while the solid line describes the model ([T]), where the homo- 
dyne efficiency 77 and the relative parasite gain 7 are kept to 
a constant value. Our model is in satisfactory agreement with 
the data, so we can be confident of that the CRB calculated 
after the matrix ( fTTj ) reliably describes the ultimate limit for 
precision. In the formulae of the Bayes rule (|6]l, we need to 
multiply several probabilities (|5]), which rapidly give a num- 
ber hardly manageable by reasonable computing power: this 
sets a limit to the number of quadrature values one can effec- 
tively use in about 800 points. In order to use larger samples, 
we have divided our data in A^;, = 10^ blocks of 200 points 
for each of the quadratures (|2]), calculated the Bayesian esti- 
mation of the discord for each block, then considered the av- 
erage weighted on the associated uncertainties. We notice that 
the a priori probabilities Q are calculated from the whole set 
of data containing Mr values: as they intervene in the eval- 
uation for each block, the overall number of resources to be 
considered is M = Ni, ■ Mt- 

The comparison between our experimental uncertainties 
and both the Cramer-Rao limit for our detection ( [T2] i and the 
quantum Cramer-Rao limit is shown in Figp] where we re- 
port the quantity iiTM — Afcr2(_D)/(i^^^)D£i(orthe analogue 
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FIG. 2: Experimental values of Gaussian quantum discord from ho- 
modyne data and Bayesian estimation. The points correspond to the 
estimated experimental values, while the solid line is the theoretical 
prediction for r;=0.62 and 7=0.73 (the value of 7 has been extracted 
from a best-fit of the points). Uncertainties are within the point size. 



quantity involving the QFI) expressed in dB. Km is the vari- 
ance of the discord estimator from homodyne data multiplied 
by the number of resources and divided by the relevant ele- 
ments of the (quantum) inverse Fisher matrix. For Km equal 
to unity we have optimal estimation. Solid points refer to 
Bayesian estimation while empty ones correspond to estima- 
tion by inversion. We notice that for low values of discord, the 
Bayesian technique provides a nearly optimal estimator for 
the chosen measurement strategy, whereas estimation by in- 
version is noisier. We also notice that the point corresponding 
to the lowest value of the discord is slightly below the quan- 
tum CRB: this confirms that for low values of the squeezing 
of the pump, the model we use is not as accurate as in other 
regimes. For increasing values, the observed variances de- 
part from the optimum by less than an order of magnitude: as 
Bayesian estimation rapidly converges to optimal, we can at- 
tribute this trend to actual variations of the value of the discord 
in the experiment, becoming more important than statistical 
fluctuations when the discord increases. 

The measurement we have adopted has the considerable ad- 
vantage of being the simplest experimental option; however, 
simplicity always comes at a price, and we do not expect it 
to deliver the best estimator for discord as established by the 
quantum CRB. In the limit of low discord, we measure a ratio 
of about 10 dB , which tells us that the price we have to pay is 
quite reasonable. The departure from the quantum CRB then 
slightly increases with discord. 

In conclusion, we have presented the experimental esti- 
mation of Gaussian quantum discord for two-mode squeezed 
state. Our scheme is based on homodyne detection assisted 
by Bayesian analysis. Our results are in good agreement with 
the theorerical model, and this allows us to perform a reliable 
precision analysis. We found that homodyne estimation shows 
about 10 dB of added noise compared to the ultimate bound 
imposed by the quantum Fisher information, with Bayesian 
analysis that slighlty improves performances for small values 




FIG. 3: The noise ratio Km as a function of discord. Circles and 
triangles correspond respectively to the quantum and the classical 
CRB. Solid points refer to Bayesian estimation, while empty ones 
correspond to estimation by inversion. The uncertainties for the in- 
version method are estimated using a MonteCarlo procedure on 10^ 
points. Notice that the number of resources for Bayesian estimation 
is M — Ni, ■ Mt, while for the inversion method A/ — Mt- 



of discord. We have also compared our results with the CRB 
for homodyne detection and found that the estimation is nearly 
optimal for small values of discord. The 

usefulness of quantum discord as a resource for quantum 
technology is a heavily debated topic, and a definitive an- 
swer may only come from experiments involving carefully 
prepared quantum states. Our results contribute to the pre- 
cise characterization of Gaussian discord and illustrate how a 
suitable data processing may decrease the uncertainty when 
optimal detection schemes are not available. 
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SUPPLEMENTAL MATERIAL 



Definition of quantum discord 



The state of a bipartite system qab is called separable if it can be produced by local operations and classical communication, 
viz. QAB — X] PkCAk ^ <^Bk, o'Ak and aBk are generic density matrices describing the states of the two subsystems. Despite the 
fact that uniquely classical information is exchanged, this procedure can neverless generate quantum correlations, as revealed by 
inspecting the mutual information of the two subsystems. This is the idea on which quantum discord is built. 

For classical variables, mutual information can be defined by the equivalent expressions 



/(A; B) = H{A) + H{B) - H{A, B), 



(14) 



and 



I{A; B) ^ H{A) - H{A\B) = H{B) - H{B\A), 



(15) 



where H{X) = — ^^ px {x) logpx (x) is the Shannon entropy of the corresponding probability distribution and the conditional 
entropy is defined as 7J(A|B) = -Y.bPB{b)J2aPA\B{a\b) \ogpA\B{a\b) = -Y.abPAB(a-,b) log pA|B(a|^), where we used 
the joint probability p^B(a, 6), the two marginal probabilities pA(a) = X^bP^sl'^' ^) ^^^ PB{b) = J2aPA-B{a,b), and the 
conditional probabilities PA|B(a|&) = p{o-, b)/p{b) and PB\A{b\a) — p{a, h)/p{a). 

The idea of quantum discord grows out of the fact that the quantum version of the mutual information of a bipartite state qab 
may be defined in two inequivalent ways. The first is obtained by the straightforward quantization of the classical expression 
([T4J, 



I{qab) = S{qa) + S{qb) - S{qab) 



(16) 



where S{q) = — Tr[g log q] is the Von-Neumann entropy of the state q and qa — '^'^b[qab]7 Qb = Tr^[p^s] are the partial 
traces over the two subsystems. On the other hand, the quantization of the expression based on conditional entropy involves the 
conditional state of a subsystem after a measurement performed on the other one. This fact has three relevant consequences: 

• the symmetry between the two subsystems is broken; 

• this quantity depends on the choice of the measurement; 

• the resulting expression is generally different from I {qab)- 



Let us denote by gAk = T^b[qab I® Pk]/pB{k) withpsik) = Tr ab[qab ^<E} Pk], the state of the system A conditioned on the 
outcome k from a measurement performed on the system B; {Pk} denotes the elements of a POVM . The quantum analogue of 
the expression ( [T5] l is then defined as the upper bound 



Ja = sup S{gA) - y^PBik)S{gAk) 
{P.} V 



(17) 



taken over all the possibile measurements. This represents the maximal Holevo information that can be achieved by using the 
subsystem A while adopting subsystem i? as a measuring device for A. 

Finally, the quantum A-discord is defined in terms of the mismatch D{qab) = I{qab) — Ja{qab)- Analogously one is led 
to define the B-discord through the entropy of conditional states of system B. 

The direct transposition of these definitions to the continuous variable realm is hindered by the complexy of the maximisation 
requested in ( [T7| ) in an infinite-dimensional space. For Gaussian states, one can introduce the notion of Gaussian quantum 
discord by restricting the maximisation to Gaussian POVMs, which lead to an analytical expression. In order to obtain an 
explicit expression, we write the covariance matrix of the state in the form 



/ a c Q\ 

a -c 

c 6 
yO -c h J 

The symplectic invariants are then given by /i=a^, l2=h^, Ij,^ — (?, I4 = det(S). We can calculate the A-discord as 

V 1 + 2V72 



(18) 



D{Y.)=h{y/ 12) - h(d_) - h{d+) + h 



(19) 



In the formula above we have introduced the symplectic eigenvalues d\ — }^ [^ ^ ^/A^— 4/4] , A=/i + / + 2 + 2/3, and h{x) 
is the binary entropy 



h{x) = (x + 1/2) log(x + 1/2) - (a; - 1/2) log(x - 1/2) 



(20) 



For symmetric states - such as those considered in the present investigation - the distinction between A-discord and B-discord 
becomes superfluous and one can simply talk about the discord of the state. 

Multiparametric quantum estimation 

Here we present the case where the estimation of more than one parameter has to be performed. We define a family of quantum 
states Qx which depends on a set of iV parameters A = {A^}, jjl = I, . . . ,N. In this case the geometry of the estimation problem 
is contained in the QFI matrix, whose elements are defined as 



H{X)^, - Tr 



9\- 



J^n^l/ + J^ljJ^L 



(21) 



and where we have introduce the Symmetric Logarithmic Derivatives (SLD) L^ corresponding to the parameter A^, as the 
selfadjoint operator that satisfies the equation 



Lfj.gx + QxL^ _ dgx 
2 ^ dX^- 

In terms of eigenvalues and eigenvectors of gx, by denoting with 9^ the partial derivative respect to A^, we have 

H{\}^,^ = 2^ -t- 2^ X 



(22) 



n^m 



an + a. 



n i^ "-m 



X i{'ipn\dfx1p„i){d^1pm\^n) + (V"™ l^z.l/'m) (^pl/'m | V'ri)) ■ 



(23) 



The QFI matrix here defined provides a lower bound (the quantum Cramer-Rao bound) on the covariance matrix 7 = (A^ A^) — 
(A^)(A^), i.e.. 



7>^H(A)- 



(24) 



In the multiparametric case this bound is not in general achievable, on the other hand, the diagonal elements of the inverse Fisher 
matrix provide achievable bounds for the variances of single parameter estimators, at fixed value of the others 

Var(A^)=7^^>^iJ(A);^i. (25) 

Let us now suppose that we are interested in the estimation of different set of parameters A = {X^, = X^{\)} which are functions 
of the previous ones. We then need to reparametrize the family of quantum states in terms of A. Since d^ = ^ Bf^^^df^ with 

B^p = dXf^/dXi, we have that 

L.^Yl ^A'-^M (26) 

and the new QFI matrix simply reads 

H = BHB^. (27) 

We consider here the case where we perform a specific indirect measurement in order to infer the values of the parameters A, 
given some measurement outcomes X — {xi, X2, . . . }. The whole measurement process can be described by the conditional 
probability p{x\X) of obtaining the value x from the measurement when the parameters have the values A. Given this object, we 
can define the Fisher information (FI) matrix whose elements are obtained as 

F.^-jdxv{^\X) |;j |;j \ (28) 

This matrix defines a bound on the covariance matrix 7 for the specific measurement we performed. In particular we are 
interested in the bound for the variance of a single-parameter, at fixed values of the others, which reads 

Var(A^) > l-^F{\)-l > j^HiX);'^ (29) 

and which in turn is always lower bounded by quantum CRB given in Eq. ( [25] l. 

Notice that if we have to reparametrize our family of states in terms of different parameters A, we can use the same formulas 

shown above for the QFI matrix, obtaining the new FI matrix as -F = BFB. 

Physical model and evaluation of quantum discord 

Here we give explicit expressions of the formulas used in the main text. A two-mode squeezed thermal state (STS) is fully 
characterized by the two parameters Ng ~ sinh s and Nt, representing, respectively, the effective amount of squeezing photons 
and thermal photons. In our experimental model, these quantities can be obtained as a function of the physical parameters 
{r, 7, 7^}, that is 



2 



where 



W 7] 2 cosh rcosh (2r7) + i?(r, 7^)^ + 277 cosh r(~2rycosh (r7) sinli r + cosh(2r7)i?(r, 77)) ^ 
Nt ~ - 1 — 1 + Y ( A(r, 7, 77) — 77 cosh^ (rj) sinh 2r) {A{r, j.rj) +77 cosh^ {r"f) sinh 2r) J (30) 

y4(r, 7, 77) = 1 — 77 + 77 cosh r cosh 2r7 + 77 sinh r (31) 

_B(r, 77) = 1 — 77 + 77 sinh r. (32) 

Notice that by varying the pump power, we change the parameter r only, while the noise parameters 7 and 77 stay constant to a 
very good level of approximation (it should depend on the mode matching only). As a result, both the effective squeezing and 
thermal photons Ng and Nt change accordingly. 
The covariance matrix of a two-mode STS can be written as 

^- - ( f: 7i ) ^33) 

\caz al2 J 



where 



a={l + 2Nt){l + 2Ns) (34) 

c^2{l + 2Nt)VNsiN, + l), (35) 



and I2 and 0-^ are respectively the 2x2 identity matrix and the PauH matrix for the z direction. Following 11261 . the quantum 
discord can be thus evaluated, obtaining 

D{N,,Nt) = 2Nt log{Nt) - 2{Nt + 1) \og{Nt + 1) - {N, +Nt + 2N,Nt) log(iV, +Nt + 2NsNt)+ 
^*^^* ^ ^^ ■ log r , .!^*,^^t^\^. .. ) +{l + Ns + Nt + 2N,Nt) log(l + N, + Nt + 2NsNt) 



l + Ns+Nt + 2NsNt ^\l + Ns+Nt + 2NsNtJ ' - w e.v .1; ^3^^ 

, Ns + 2N,Nt + il + Nt)\ fN, + 2NsNt + (1 + Ntf 

log 



1 + A^, + iVi + 27V,iVt °V l + Ns+Nt + 2N,Nt 



expressed as a function of the effective parameters Ng and Nf. By using Eqs. (30i, one can easily obtain the discord as a 
function of the physical parameters r, 7 and 77. 



Monte Carlo evaluation of uncertainties. 

In our experiment, Mq measurement outcomes are recorded for each one of the four quadratures. The quadratures 
Q^^) and Q^*) form a set of 2Mq squeezed quadratures measurements, as well as Q^^^ and Q^''^ form a set of 2Mq anti- 
squeezed quadratures measurements. From those two sets of experimental data, we compute the variances cr^((5sq) and 
c^(Qasq)- Assuming that those estimated variances follow a Gaussian distribution, the variance of their estimation is given 
by Vai- (a2(Q,q/,3q)) ^2a\Q,^/,,^)/[2Mq). 

The expressions for a^{Qsq) and (T^(Qasq) are then inverted to obtain Ns and Nt as function of (T^((3sq) and ^^(Qasq)- 
The mean value and the associated uncertainties are determined by a Monte Carlo simulation of 10^ experiments. For each 
experiment, the values iT^((5sq) and o'^(Qasq) of the squeezed and anti-squeezed variances are randomly chosen from two 
Gaussian distributions respectively of mean values a'^{Qsq) and (T^{Qasq), and variances Var (o'^(Qsq)) and Var (^^(Qasq))- 
The values of Ng and Ni are then computed using those random values. The experimental estimate A^"^ and iV^"^ are finally 
obtained by taking the mean of the 10^ values of Ng and Nt, whereas their uncertainties (y'^{Nf") and (j'^{Nf") are obtained 
by computing the variance of the 10^ values of Ng and Nf 

These values can be used in the expression for discord ([36| to calculate its value 13'"^, and the uncertainty (j^{D™^) by using 
a similar MonteCarlo method. 



